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In this article, we propose a simple rule that generates scale-free networks with very large cluster- 
ing coefficient and very small average distance. These networks are called Random Apollonian 
Networks(RAN) as they can be considered as a variation of Apollonian networks. We obtain the 
analytic results of power-law exponent 7 = 3 and clustering coefficient C = ^ — 361n| ~ 0.74, 
which agree very well with the simulation results. We prove that the increasing tendency of average 
distance of RAN is a little slower than the logarithm of the number of nodes in RAN. Since most real- 
life networks are both scale-free and small-world networks, RAN may perform well in mimicking the 
reality. The RAN possess hierarchical structure as C{k) ^ k^^ that in accord with the observations 
of many real-life networks. In addition, we prove that RAN are maximal planar networks, which are 
of particular practicability for layout of printed circuits and so on. The percolation and epidemic 
spreading process are also studies and the comparison between RAN and Barabasi- Albert (B A) as 
well as Newman- Watts(NW) networks are shown. We find that, when the network order A''(the 
total number of nodes) is relatively small(as A*" ~ 10*), the performance of RAN under intentional 
attack is not sensitive to A'^, while that of BA networks is much affected by A'^. And the diseases 
spread slower in RAN than BA networks during the outbreaks, indicating that the large clustering 
coefficient may slower the spreading velocity especially in the outbreaks. 

PACS numbers: 89.75.Hc, 64.60.Ak, 84.35.-l-i, 05.40.-a, 05.50-fq, 87.18.Sn 



I. INTRODUCTION 

Many social, biological, and communication systems 
can be properly described as complex networks with 
nodes representing individuals or organizations and edg^ 

mimicking the interactions among themiE mm nn 

0, i, H El m m m EI. Examples are numerous: 
these include the Internetjll El [13, the World Wide 
WebEi,Ei|2i|2l|22, social networks of acquaintance 
or other relations between individuals 123. 125l l2^ 
IIE I30I I31L Is^ ■ metabolic networks |34 l35l ISa . 
l37l l38l , food webs 39, MElElJi^^ many 
others 149, 

ubiquity of complex networks inspires scientists to con- 
struct a general model. In the past 200 years, the study of 
topological structures of the networks used to model the 
interconnection systems has gone through three stages. 
For over a century, there is an implicit assumption that 
the interaction patterns among the individuals can be 
embedded onto a regular structure such as Euclidean 
lattices, hypercube networks, and so onjsM Isol l60l l6ll|. 
Since late 1950s mathematicians began to use random 
graphs to describe the interconnections, this is the sec- 
ond stage m m 111 m S mill m . in the past few 
years, with the computerization of data acquisition pro- 
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cess and the availability of high computing powers, scien- 
tists have found that most real-life networks are neither 
completely regular nor completely random. The results 
of many experiments and statistical analyses indicate 
that the networks in various fields have some common 
characteristics, the most important of which are called 
small-world effect. 70. 71. 72. J4, and scale-free 
propertylzl III EflHr^ 

In a network, the distance between two nodes is defined 
as the number of edges along the shortest path connect- 
ing them. The average distance L of the network, then, 
is defined as the mean distance between two nodes, av- 
eraged over all pairs of nodes. The average distance is 
one of the most important parameters to measure the ef- 
ficiency of communication networks. For instance, in a 
store-forward computer network, probably the most use- 
ful measurement of its performance is the transmission 
delay (or time delay) encountered by a message travelling 
through the network from its source to destination; this 
is approximately proportional to the number of edges a 
message must pass through. Thus the average distance 
plays a significant role in measuring the transmission de- 
lay. The number of the edges incident from a node x 
is called the degree of x, denoted by k{x). Obviously, 
through the k{x) edges, there are k{x) nodes that are 
correlated with x; these are called the neighbor-set of x, 
and denoted by A{x). The clustering coefficient C{x) of 
node X is the ratio between the number of edges among 
A{x) and the total possible number, the clustering co- 
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FIG. 1; An Apollonian packing of disks within a circle fTi^ . 



efficient C of the whole network is the average of G{x) 
over aU x. Experiments indicate that most real-life net- 
works have much smaller average distance (as L r^hiN 
where N is the number of nodes in the network) than 
the completely regular networks and have much greater 
clustering coefficient than those of the completely ran- 
dom networks. Therefore they should not be treated 
as either completely regular or random networks. The 
recognition of small-world effect involves the two factors 
mentioned above: a network is called a small-world net- 
work as long as it has small average distance and great 
clustering coefficient. Another important characteristic 
in real-life networks is the power-law degree distribution, 
that is p{k) oc , where k is the degree and p{k) is the 
probability density function for the degree distribution. 
7 is called the power-law exponent, and usually between 
2 and 3 in real-Hfe networks^ H Q. This power-law 
distribution falls off much more gradually than an expo- 
nential one of a completely random network, allowing for 
a few nodes of very large degree to exist. Networks with 
power-law degree distribution are referred to as scale- free 
networks, although one can and usually does have scales 
present in other network properties. 

From 1998 much attention has been focused on how 
to model complex network. One of the most well-known 
models is Watts and Strogatz's small- world network (WS 
network), which can be constructed by starting with a 
regular network and randomly moving one endpoint of 
each edge with probability p 71] . Another popular model 
was proposed independently by Monasson|7^ and by 
Newman and Watts (NW Networks) 'tS] , where no edges 
are rewired. Some variations of the small-world model 
have been proposed. Several authors have studied the 
model in dimension higher than one and obtained simi- 
lar results to the one dimension case|73, 113, ^1, 82, 83]. 
Another kind of models in which shortcuts preferentially 
join nodes that are close together on the underlying lat- 



tice also have been studied '84','85','8^. Very recently, Zhu 
et. al. proposed a so-called directed dynamical small- 
world model, in which the network structure is affected 
by the processes upon the networkfTsj. 

Another significant model is Barabasi and Albert's 
scale- free network model (BA network) [t^ E3 ' which is 
very much similar to Price'sfT^IX^]- The BA model sug- 
gests that two main ingredients of self-organization of a 
network in a scale-free structure are growth and prefer- 
ential attachment. These point to the facts that most 
networks grow continuously by adding new nodes, which 
are preferentially attached to existing nodes with large 
number of neighbors. The subsequent researches on var- 
ious processes taking place upon complex network^j_such 
as percolation 8^181 El IM I M lM[ 9d |9 ^^ 96, 97 , 

ai, epidemic processes lggniMTim^l^^ , 1104 il05 , 
lOa . cascade processes [iOtI llOSl Il09l iTTa iTTl iTTa 1113 . 

and so on, indicate that the scale-free degree 
distribution plays the most crucial role rather than small- 
world effect. Therefore, in the recent two or three years, 
the study of modelling complex networks focuses on re- 
vealing the underlying mechanism of power-law degree 
distribution. Roughly, these models for scale-free net- 
works can be classified into 3 main scenarios. The first 
one is related to the models of human behavior |ll6llll'^ 
and was introduced in a network version under the name 
"pre ferential attachment" mentioned abo yel?^ 17^ IllSl 

ITi9llTM IiHIT2aiMITMITMIiMIT27ili2 ^^ 

The second class of models is where a scale-free dis- 
tribution appears as a result of a balance between a 
modelled tendency to form hubs against an entropic 
pressure towards a random network with an exponen- 
tial degree distribution 131, 132, 133]. The third one is 
the self-organized models that lead to p ower-law degree 
distribution[l3l [HI [isi [isl [Hi [iSl [m^ . 

In the recent several months, a few authors have 
demonstrated the use of pure mathematical objects and 
methods to construct scale-free networks. The first in- 
teresting instance is the so-called integer networks jl41j, 
in which the nodes represent integers and two nodes x 
and y are linked by an edge if and only if x is divided ex- 
actly by ?/ or y is divided exactly by x, where x and y are 
nonzero integers. Zhou and Wang et al have studied the 
statistical properties of these networks and demonstrated 
that they are scale- free networks of large clustering coeffi- 
cient. Another significant instance is Apollonian Net- 
virorks(AN) introduced by Andrade et al 142]. In our 
opinion, Apollonian networks may be not the networks 
of best performance, but assuredly the most beautiful 
ones we have ever seen. Another related work is owed to 
Dorogovtsev and Mendes et. al., in which the determin- 
istic networks, named pseudofract als, are o btained by 
random attachment aiming at edges [T43.ll44j . 

In this article, we propose a simple rule that generates 
scale-free networks with very large clustering coefficient 
and very small average distance. These networks are 
called Random Apollonian Networks (RAN), since 
they can be considered as a variation of Apollonian 
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FIG. 2: The development of the 2D ApoUonian network inside 
the interstice between three mutually touching disks, as the 
number of generations increase. In e a.ch p icture, the network 
is overlaid on the underlying packing |l48j| . 

networks |l45j . We discuss the difference between RAN 
and AN in detail(in addition, the difference between RAN 
and BA networks), and demonstrates that RAN perform 
much better than BA networks in some aspects. 

This article is organized as follow: In section 2, the 
Apollonian networks is introduced, including its fasci- 
nating properties and its shortcomings. In section 3, the 
rule that generates RAN is described in detail. In sec- 
tion 4, we give both the simulation and analytic results 
about the statistical characteristics of RAN, including 
the scale- free property and the small- world effect, where 
the detailed analytic processes are shown in appendices. 
In section 5, we prove that RAN are the maximal planar 
networks. In section 6, The percolation and epidemic 
spreading process are also studies and the comparison 
between RAN and BA networks as well as NW networks 
are shown. Finally, we sum up this article and discuss 
the relevance of RAN to the real world in section 7. 



exponent of about 1.3 |l47j |. Of course, these added disks 
cannot fill all of the space in the interstices, but instead 
give rise to three smaller interstices. In the second gener- 
ation, further disks are added inside all of these new in- 
terstices, which again touch the surrounding disks. This 
process is then repeated for successive generations. If we 
denote the number of generations by t, where t = cor- 
responds to the initial configuration, as t — > cxd the space- 
filling Apollonian packing is obtained as shown in figure 
1. Apollonian packing can be used as a basis of a net- 
work, where each disk is a node in the network and nodes 
are connected if the corresponding disks are in contact. 
We call this contact an "Apollonian Network" . Figure 
2 shows how the network evolves with the addition of 
new nodes at each generation. For each new disk added, 
three new interstices in the packing are created, that will 
be filled in the next generation. Equivalently, for each 
new node added, three new triangles are created in the 
network, into which nodes will be inserted in the next 
generation. 

Doye et. al. have studi ed the properties of Apollonian 
networks detailedly JMS] , and shown the degree distri- 
bution p{k) oc k~'^ , average length I oc (InN)^ , where 
7 = 1 + 1^ « 2.585, (3 « 0.75 and N is the order [l49t 
of the network, in other words, Apollonian networks are 
scale-free and display small-world effect. It is worth re- 
marking that the clustering coefficient C is close to 0.828, 
much larger than that of BA networks, in the limit of 
large TV. 

Andrade et. al. have also found many peculiar results 
about some well-known models upon Apollonian net- 
works, including percolation, electrical conduction, and 
magnetic model (142. .14(| . 



III. RANDOM APOLLONIAN NETWORKS 



II. BRIEF INTRODUCTION OF APOLLONIAN 
NETWORKS 

Ap ollonian networks, introduced by Andrade et al 
[142'| , are derived from the problem of space- filling pack- 
ing of spheres according to th e ancient Greek mathemati- 
cian ApoUonius of Perga |l47j . To produce an Apollonian 
packing, we start with an initial array of touching disks, 
the interstices of which are curvilinear triangles. In the 
first generation disks are added inside each interstice in 
the initial configuration, such that these disks touch each 
of the disks bounding the curvilinear triangles. The po- 
sitions and radii of these disks can easily be calculated, 
and the circle size distribution follows a power-law with 



Random Apollonian network starts with a triangle con- 
taining three nodes marked as 1, 2 and 3. Then, at each 
time step, a triangle is randomly selected, and a new 
node is added inside the triangle and linked to the three 
vertices of this triangle. The sketch maps for the network 
growing process are shown in figure 3. 

It is clear that, at time step t, our network is of order 
N — t + 3. Using this simple rule, one can get ran- 
dom Apollonian network of arbitrary order as he like. 
Note that the randomicity is involved in our model(that 
is why we call these networks random Apollonian net- 
works), the analytic approaches are completely different 
from the earlier studies on Apollonian networks. 
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FIG. 3: The sketch maps for the network growing process. 
The four figures show a possible growing process for RAN at 
time t = 1(a), t = 2(b), t = 3(c) and t = 4(d). At time step 1, 
the 4th node is added to the network and linked to node 1, 2 
and 3. Then, at time step 2, the triangle A134 is selected, the 
5th node is added inside this triangle and linked to node 1, 3 
and 4. After that, the triangles A234 and A 124 are selected 
at time step 3 and 4, respectively. And node 6 and 7 are 
added inside these two triangles respectively. The figure 3d 
shows a random Apollonian network of order 7. Keep on the 
similar iterations, one can get RAN of any orders as he like. 



FIG. 4: Degree distribution of RAN, with N = 640000(black 
squares), N = 320000(red circles) and iV = 160000(blue tri- 
angles). In this figure, P{k) denotes the number of nodes of 
degree k. The power-law exponent 7 of the three probability 
density function are 7640000 = 2.94±0.04, 7320000 = 2.92±0.05 
and 7160000 = 2.92 ± 0.06, respectively. The average expo- 
nent of them are 2.93. The inset shows Degree distribution of 
RAN, with iV = 80000(black squares), N = 40000(red circles) 
and A'' = 20000(blue triangles). The exponents are 7540000 = 
2.91 ± 0.07, 7320000 = 2.90 ± 0.07 and 7160000 = 2.90 ± 0.09. 
The mean value is 2.90. The two dash lines have slope -3.0 
for comparison. 



IV. STATISTICAL CHARACTERISTICS OF 
RANDOM APOLLONIAN NETWORKS 

A. The scale-free property 

As we have mentioned above, the degree distribution 
is one of the most important statistical characteristics 
of networks. Since a majority of real-hfe networks are 
scale-free networks, whether the networks are of power- 
law degree distribution is a criterion to judge the vahdity 
of the model. In this subsection, we will give the simula- 
tion and analytic results on random Apollonian networks' 
degree distribution. 

Note that, after a new node is added to t he networks, 
the number of triangles increases by 2 jl5n| |. Therefore, 
we can immediately get that when the networks are of 
order TV, the number of triangles are: 

A^A = 2(7V-3) + 1 = 2iV- 5 (1) 

Let N^^ denote the number of triangles containing the 
ith node, the probability that a newly added node will 
link to the ith node is N'^/N/\. Apparently, except the 
node 1, 2 and 3, A''^ is equal to the degree of the ith 
node: N\ = ki. Therefore, we can write down a rate 
eauation |ll8j for the degree distribution. Let n{N, k) be 
the number of nodes with degree k when N nodes are 



present, now we add a new node to the network, n{N, k) 
evolves according to the following equation: 

n{N+l, fc+1) = n{N, k)^+n{N, fc+l)(l-^) (2) 

When N is sufficient large, n{N^ k) can be approximated 
as Np{k), where p{k) is the probability density function 
for the degree distribution. In terms of p{k), the above 
equation can be rewritten as: 

iVA iVA 

(3) 

Using Equ.(l) and the expression p{k -f- 1) — p{k) — j^, 
we can get the continuous form of Equ.(3): 

, dp 3N — 5 ,, , „ 

This lead to p{k) oc k^f with 7 = (3A^ - 5)/A^ w 3 for 
large N. 

In figure 4, we report the degree distribution for 
N = 640000,320000,160000,80000,40000 and 20000. 
The simulation results agree very well with the analytic 
one. 

By the way, some readers may think the RAN are al- 
most the same as BA networks. Indeed, the two ingredi- 
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FIG. 5: The clustering coefficient of RAN (black squares) and 
BA networks(red circles). In this figure, one can find that 
the clustering coefficient of RAN is almost a constant a little 
smaller than 0.77, which accord with the analytic very well. 
The dash line represents the analytic result 0.74. It is clear 
that, the clustering coefficient of BA networks is much smaller 
than that of RAN. The inset shows that the clustering coef- 
ficient of BA networks decreases with the increasing of net- 
work order, following approximately a power law C ~ A'^""'^'', 
which is quite different form the real-life networks. All the 
data are obtained by 10 independent simulations. 



ents "growth" and "preferential attachment" are in com- 
mon and they have almost the same power-law exponent, 
but the further simulation and analytic results will show 
that RAN and BA networks are essentially different in 
some aspects and RAN may be closer to reality rather 
than BA networks. 



B. The small- world effect 

The small- world effect consists of two properties: large 
clustering coefficient and small average distance. In this 
subsection, we will give both the simulation and analytic 
results about the two properties and prove that RAN are 
small- world networks. 

At first, let us calculate the clustering coefficient of 
RAN. As we mentioned in section 1, for an arbitrary 
node X, the clustering coefficient C{x) is: 



C{x) 



2E{x) 



k{x){k{x) ~ 1) 



(5) 



where E{x) is the number of edges among node a:'s 
neighbor-set A{x), and k{x) = \A{x)\ is the degree of 
node X. The clustering coefficient C of the whole net- 
work is defined as the average of C{x) over all nodes. 

By means of theoretic calculation(see Appendix A 
for details), we obtain the clustering coefficient of RAN 



FIG. 6: The dependence between the average distance L 
and the order of RAN. One can see that L increases very 
slowly as TV increases. The inset exhibits the curve where L 
is considered as a function of In A, which is well fitted by 
a straight line. The curve is above the fitting line when A'^ 
is small(1000 < N < 3000) and under the fitting line when 
N is large(A'^ > 4000), which indicates that the increasing 
tendency of L can be approximated as InA and in fact a little 
slower than InA. All the data are obtained by 10 independent 
simulations. 



with large order N as: 



C^-- 36ln- 



0.74 



(6) 



Figure 5 shows the simulation results about the clus- 
tering coefficient of RAN, which agree very well with the 
analytic one. It is remarkable that, the clustering coef- 
ficient of BA networks is very small and decreases with 
the increasing of network order, following approximately 
a power law C ^ N"^'"^^ our simulation, it is found 
to be about 0.73). The simulation about clustering coef- 
ficient of BA network can also be seen in figure 5. Since 
the data-flow patterns show a large amount of cluster- 
ing in interconnection networks, the RAN may perform 
better than BA networks. In addition, the demonstration 
exhibits that most real-life networks have large clustering 
coefficient no matter how many nodes they have. That 
is agree with the case of RAN but confiict with the case 
of BA networks, thus RAN may be more appropriate to 
mimic the reality. 

In addition, many real networks including Internet, 
World Wide Web, and the actor network, are character- 
ized by the existence of hierarchical structure 151, 15^ 
Il53j| . which can usually be detected by the negative cor- 
relation between the clustering coefficient and the de- 
gree. The BA networks , which does not possess hierar- 
chical structure, is known to have the clustering coeffi- 
cient C{x) of node x independent of its degree fcfa:) |l5lj . 
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FIG. 7: Random site-percolation transition for RAN(blue 
triangles), NW(red circles) and BA(black squares) networks. 
Plotted is the fraction of nodes that remain in the giant com- 
ponent after random breakdown of a fraction p of all nodes, 
riijp), as a function of p. All the networks using for numerical 
study are of order A'^ = 10000 and (fc) = 6. Clearly, the scale- 
free networks(BA and RAN) are much more resilient than 
networks of single-scale(NW), which agrees with the well- 
known conclusion 87, 90j. When p < 0.3, the performances 
of RAN and BA networks under random failures are almost 
the same; when 0.3 < p < 0.6, RAN are little resilient than 
BA networks; when p becomes even larger, BA networks get 
obviously more resilient than RAN. The critical thresholds of 

0.85 and pf ^ 0.95, 
All 



RAN and BA networks are 

which will approaches 1 as the networks grow in size 90] . 
the data are averaged over 100 independent simulations. 



while the RAN has been shown to have C(fc) ~ fc^^(see 
Eq.(A2), in ac cord with the observations of many real-life 
networks |l5lj .) 

In succession, let's discuss the average distance of 
RAN. By means of theoretic approximate calculation(see 
Appendix B for details), we prove that the increasing 
tendency of L{N) is a little slower than InA^. In figure 
6, we report the simulation results on average distance of 
RAN, which agree with the analytic result well. 

In respect that the random Apollonian networks are 
of very large clustering coefhcient and very small average 
distance, they are not only the scale- free networks, but 
also small- world networks. Since most real-life networks 
are both scale- free and small- world networks, RAN may 
perform better in mimicking reality rather than WS net- 
works and BA networks. 



V. THE PLANARITY OF RANDOM 
APOLLONIAN NETWORKS 

There are many practical situations in which it is im- 
portant to decide whether a given network is planar, and 



if so, to then find a planar embedding of the network. 
For example, a very large scale integrated (VLSI) designer 
has to place the cells on printed circuit boards according 
to several designing requirements. One of these require- 
ments is to avoid crossings since they may lead to unde- 
sirable signals. One is, therefore, interested in knowing if 
a given electrical network is planar, where the nodes cor- 
respond to electrical cells and the edges to the conductor 
wires connecting the cells. 

A network is a planar network if it can be drawn in the 
plane in such a way that no two edges intersect. Putting 
it a little more rigorously, it is possible to represent it 
by a drawing in the plane in which the nodes correspond 
to distinct points and the edges to simple Jordan curves 
connecting the points of its end-points. In this drawing 
every two curves are either disjoint or meet only at a 
common end-point. The above representation of a graph 
is said to be a plane networks. 

In some places, the networks will perform better when 
they have more edges. Therefore, how to add more edges 
into a network but keeping it a planar network is a prac- 
tical and interesting problem. According to the rule that 
generates RAN, one can immediately find that RAN are 
planar networks. Hereinafter, we will show that the RAN 
are maximal planar networks, which are the planar net- 
works with fixed order who has maximum edges. 

If we omit the nodes and edges of a planar networks 
from the plane, the remainder falls into connected com- 
ponents, called faces [bJ. Clearly, each plane network 
has exactly one unbounded face and each edge is in the 
boundary of two faces. If we draw the graph of a convex 
polyhedron in the plane, then the faces of the polyhedron 
clearly correspond to the faces of the plane networks. 
This lead to the Euler's polyhedron theorem or simply 
Euler's formula j59ll6Cll 16 1| that if a connected plane net- 
works has n nodes, m edges and / faces, then: 



n — m + / = 2 



(7) 



Furthermore, denote fi the number of faces having ex- 
actly i edges in their boundaries. Clearly, 



(8) 



And since each edge is in the boundary of two faces, we 
have: 



^ if^ = 2m 



(9) 



Combine Equ.(8) and Equ.(9) and note that each face 
has at least 3 edges, we have: 

2m>3^/, = 3/ (10) 

i 

Then, using Euler's formula, one can obtain that: 



TO < 3n — 6 



(11) 
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That is to say the maximum number of edges for a planar 
network with order n is 3rt — 6. Apparently, the random 
Apollonian network with order N has 3iV — 6 edges on the 
beam, thus all the RAN are maximal planar networks. 



VI. PERCOLATION AND EPIDEMIC 
SPREADING ON RANDOM APOLLONIAN 
NETWORKS 

As we mentioned above, close to many real-life net- 
works, random Apollonian networks are both scale-free 
and small-world. Therefore, it is worthwhile to investi- 
gate the processes taking place upon RAN and directly 
compare these results with just small- world networks(like 
WS and NW networks) and just scale-free networks(like 
BA networks). These comparisons may give us deep sight 
into dynamic properties of networks. 

In this section, we will exhibit some simulation results 
on two deeply studied ones, percolation and epidemic 
spreading process. Since the WS networks may be un- 
connected, wc will use NW networks as exemplifications 
of small- world networks. 



A. Percolation 

In the year 2000, Albert et. al. have raised the 
questions of random failures and intentional attack on 
networksll^. Part of these questions can be equally 
considere d as site -percolation or bond-percolation on 
networks 154, 155]. In this subsection, upon RAN, BA 
and NW networks, we will investigate two kinds of site- 
percolation, random site-percolation(RSP) and pref- 
erential site-percolation (PSP), which correspond to 
random failures and intentional attack aiming at nodes, 
respectively. 

When such networks are subject to random 
breakdowns-a fraction p of the nodes and their in- 
cident edges are removed randomly-their integrity might 
be compromised: when p exceeds a certain threshold, 
p > Pc, the network disintegrates into smaller, discon- 
nected fragments. Below that critical threshold, there 
still exists a connected cluster named giant component 
that spans the entire system(its size is proportional 
to that of the entire system). In figure 7, we report 
the simulation results for RSP upon RAN, BA and 
NW networks. RAN and BA networks are obviously 
more resilient than NW networks under random failures 
that agrees with the well-known conclusion [s^ [oS. 
According to the results obtained by Cohen et. al.|90l|. 
RAN and BA networks do not have nontrivial critical 
threshold pc < 1 in the limit of large network order 
N oo. However, for finite network order and large p, 
RAN are frailer than BA networks. 

Figure 8 shows the performances of RAN, BA and NW 
networks under intentional attack, which means the re- 
moval of nodes is not random, but rather nodes with 
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FIG. 8: Preferential site-percolation for RAN(red circles), 
BA(black squares) and NW(green triangles) networks(from 
left to right) with = 10000 and (A;) = 6 fixed. Plotted 
is the fraction of nodes that remain in the giant component 
under intentional attack of a fraction p of all nodes, ^(p), as 
a function of p. The critical threshold of RAN, BA and NW 
networks are p^"^^ « 0.03, ^ ^ 0.3, and p^^ « 0.65, 
respectively. All the data are averaged over 100 independent 
simulations. 



the greatest degree are targeted first. One can find that 
the scale-free networks are far more sensitive to sabotage 
of a small fraction of the nodes, leading support to the 
view of Albert and Cohen et. al.jUEl- Although we 
know the critical threshold for PSP upon scale-free net- 
works will decay to zero as the increasing of network order 
as limjv^ooPc = 0, the notable difference between RAN 
and BA networks when N is relatively small surprises 
us. For very large N{a,s N ^ 10^ or even larger), the 
performances of RAN and BA are almost the same with 
Pc < 0.03 joj. However, the susceptivity to order of RAN 
and BA networks are completely different. BA networks 
are very sensitive to network order, for N — 10000, its 
critical threshold is about ten times than the asymptotic 
value. And RAN are almost impervious to the change of 
network order. Since many real-life networks are of order 
in the range 10'^ to 10^, this finding may be valuable in 
practicability. 

Why RAN and BA networks display completely differ- 
ent susceptivity to order changing, does it owe to the 
difference of clustering structure or hierarchical struc- 
ture? This problem puzzles us much. To study the 
finite-size effect for PSP upon scale-free networks in de- 
tail and to use network model with tunable clustering 
coefficient 107] may reveal some news, which will be con- 
sidered in the future. 
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B. Epidemic Spreading Process 

Recent studies on epidemic spreading in complex net- 
works indicate a particular relevance in the case of net- 
works characterized by various topologies that in many 
cases present us with new epidemic propagation sce- 
narios such as the absence of any epidemic thresh- 
old below w hich the infection cannot initiate a major 
outbreak [9^ llOOj . The new scenarios are of practical in- 
terest in computer virus diffusion and the spreading of 
diseases in heterogeneous populations. However, most 
previous studies have been focused on the stationary 
properties of endemic states or the final prevalence(i. e. 
the number of infected individuals) of epidemics. For the 
sake of protecting networks and finding optimal strate- 
gies for the deployment of immunization resources, it is 
of practical importance to study the dynamical evolution 
of the outbreaks, which has been far less investigated be- 
fore. Barthelemy et. al. reported both the analytic and 
numerical results of velocity of epidemic outbreaks in B A 
networks, which leaves us very sh ort r esponse time in 
the deployment of control measures |l05j . We have stud- 
ied the same process in weighted scale-free networks and 
demonstrated that the larger dispersion of weight of net- 
works resul ts in slower spreading, which may be a good 
news for us |l06j| . In this subsection, we intend to study 
how the connectivity pattern(i. e. topological structure) 
affects the epidemic spreading process in the outbreaks. 
Numerical simulations about BA, NW and RAN net- 
works are drawn, which may give us more comprehensive 
sight into the corresponding dynamic behavior. 

In order to study the dynamical evolution of epidemic 
outbreaks, we shall focus on the susceptible-infected(SI) 
model in which individuals c an be in tw o discrete states, 
either susceptible or infected ll56l Il57l |. Each individual 
is represented by a node of the network and the edges 
are the connections between individuals along which the 
infection may spread. The total population(the network 
order) N is assumed to be constant thus if S{t) and I{t) 
are the number of susceptible and infected individuals at 
time i, respectively, then 

N^S{t) + I{t) (12) 

In the SI model, the infection transmission is defined by 
the spreading rate A at which each susceptible individual 
acquires the infection from an infected neighbor during 
one time step. In this model, infected individuals arc as- 
sumed to be always infective, which is an approximation 
that is useful to describe early epidemic stages in which 
no control measures are deployed. According to Eq.(12), 
one can easily obtain the probability that an susceptible 
individual x will be infected at the present time step is: 

A^(t) = 1-A^(^''-^) (13) 

where 9{x,t — 1) denotes the number of infected individ- 
uals at time step t — 1 in. x's neighbor-set A{x). 

We start by selecting one node randomly and assume it 
is infected. The diseases or computer virus will spread in 




200 400 



t (time steps) 



FIG. 9: Average density of infected individuals versus time in 
networks with order TV — 10000 and average degree (k) — 6 
fixed. The black, blue and red curves(from bottom to top) 
correspond to the case of NW, RAN and BA networks re- 
spectively. The NW networks are of z = 1 and 4> = 4 x 10~*, 
thus (k) ^ 2z + <j)N — 6 .158| (see also the accurate definitions 
of z and (f> in reference |82|). The spreading rate is A = 0.01. 
All the data are averaged over 10^ independent runs. 




200 400 600 800 1000 



t (time step) 



FIG. 10: The spreading velocity vs time. The network order 
and spreading rate are the same with figure 9. The blue, red 
and black curves correspond to the case of NW, RAN and BA 
networks respectively. The spreading velocity reaches a peak 
quickly. Before the peak-time, the spreading velocities of the 
three kinds of networks satisfy the inequality > > 
vf^ . All the data are averaged over 10^ independent runs. 
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the networks in according with the rule of Eq.(13). In fig- 
ure 9, we plot the average density over 1000 independent 
runs in RAN, BA and NW networks with TV = 10000 and 
(k) = 6 fixed. Obviously, all the individuals will be in- 
fected in the limit of long time as l±mt—,oci{t) = 1, where 
*(^) = "TT denotes the density of infected individuals at 
time step t. More over, the simulation results indicate 
that the diseases spread more quickly in BA networks 
than RAN, as well as in RAN than NW networks. To 
make the outcome more clear, we have calculated the 
diseases spreading velocity, which is defined as: 



di{t) I{t)-I{t-1) 



dt 



N 



(14) 



Figure 10 shows the spreading velocity vs time in RAN, 
BA and NW networks with N = 10000 and (fc) = 6. The 
spreading velocity reaches a peak quickly. Before the 
peak-time, the spreading velocities of the three kinds of 
networks satisfy the inequality: 



BA RAN NW 



(15) 



The result that diseases spread more quickly in RAN 
and BA networks than in NW networks is easy to be 
understood as the well-known conclusion: boarder de- 
gree distri butio n will speed up the epidemic spreading 

process [ggllTooj. 

Why the diseases spread more quickly in BA networks 
than RAN is a very interesting question. We argue that 
the larger clustering coefficient may slow down the epi- 
demic spreading process especially in the outbreaks. For 
an arbitrary edge e, containing two nodes x and y, ob- 
viously, the distance between x and y is d{x,y) — 1. 
Remove the edge e from quondam network, then, the 
distance between x and y will increase d'{x, y) > l(if the 
removal of e makes x and y disconnected, then we set 
d'{x,y) = N). The quantity d'{x,y) can be considered 
as edge e's score s(e) = d'{x, y) > 2, denoting the number 
of edges the diseases must pass through from x to y or 
form ?/ to x if they do not pass across e. If s(e) is small, 
then e only plays a local role in the epidemic spreading 
process, else when s(e) is large, e is of global importance. 
For each edge e, if it does some contribution to clustering 
coefficient, it must be contained in at least one triangle 
and s(e) = 2. Therefore, networks of larger clustering co- 
efficient have more local edges. RAN and BA networks 
are two extreme ones of scale-free networks. In RAN, all 
the edges are of score 2; while in BA networks, almost all 
the edges are of score larger than 2 because the clustering 
coefficient of BA networks will decay to zero quickly as 
A'' increases. Consequently, diseases spread more quickly 
in BA networks than RAN. 

The above explanation is qualitative and rough, to 
study the process upon networks with tunable clustering 
coefficient "iQTj may be useful for the present problem, 
which will be one of the future works. 



VII. CONCLUSIONS 



In conclusion, in respect that the random Apollonian 
networks are of very large clustering coefficient and very 
small average distance, they are not only the scale-free 
networks, but also small- world networks. Since most 
real-life networks are both scale-free and small-world net- 
works, RAN may perform better in mimicking reality 
rather than WS networks and BA networks. In addi- 
tion, RAN possess hierarchical structure that in accord 
with the observations of many real networks and we pro- 
pose an analytic approach to calculate clustering coef- 
ficient. Since in the earlier studies, only few analytic 
results about clusteri ng coeffic i ent of net works with ran- 
domicity are report ed |l5llll59lll60lll6lj . we believe that 
our work may enlightened readers on this subject. 

Further more, we briefly introduce the conception of 
planar network(it is also called "planar graph" in math- 
ematical language), and prove that RAN are maximal 
planar networks, which are of particular practicability for 
layout of printed circuits and so on. Although whether a 
network is planar or not is a natural and important ques- 
tion that attracts much attention for mathematicians, it 
seems not interesting for physicists and almost no perti- 
nent results are reported in the earlier studies on complex 
networks. But in fact, many real-life networks are planar 
networks by reason of technical or natural requirements, 
such as layout of printed circuits, river networks upon 
the earth's surface, vas networks clinging to cutis, and 
so forth. Since the planar networks have some graceful 
characteristics that can not be found in non-planar ones, 
researchers ought to pay more attention to networks' pla- 
narity. We wish our abecedarian work to stimulate physi- 
cists thinking more of planarity. 

The percolation and epidemic spreading process are 
also studies and the comparison between RAN and BA 
as well as NW networks are shown. In percolation model, 
we find that, when the network order TV is relatively 
small (as N ~ 10*), the performance of RAN under in- 
tentional attack is not sensitive to N, while that of BA 
networks is much affected by A'^. In epidemic spread- 
ing process, the diseases spread slower in RAN than BA 
networks during the outbreaks, indicating that the large 
clustering coefficient may slower the spreading velocity 
especially in the outbreaks. We give some qualitative ex- 
planation about how the clustering structure affect the 
spreading process in the outbreaks, but Why RAN and 
BA networks display completely different susceptivity to 
order changing is even a problem puzzling us. Those 
simulation results suggest that the clustering structure 
may affect the dynamical behavior upon networks much. 
Since many real-life networks are of great clustering coef- 
ficient, to study the proc ess upon networks with tunable 
clustering coefficient jlOTf is significant. 
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APPENDIX A: MORE DETAILS FOR 
CALCULATION OF CLUSTERING 
COEFFICIENT 

At first, let us consider the clustering coefficient C{x) 
of an arbitrary node x except the nodes 1, 2 and 3. At 
the very time when node x is added to the network, it 
is of degree 3 and E{x) = 3. After, if the degree of 
node X increases by onc(i.e. a new node is added to be a 
neighbor of x) at some time step, then E{x) will increases 
by two since the newly added node will link to two of the 
neighbors of node x. Therefore, we can write down the 
expression of E{x) in terms of k{x): 



E{x) = 3 + 2{k{x) - 3) = 2k{x) - 3 



(Al) 



Using Equ.(5), we can get the clustering coefficient of 
node x as: 



C{x) 



2(2fc(x) - 3) 
k{x){k{x) - 1) 



(A2) 



2 ^ 2,k- 3 2^3 1 
^=NY.k,{k,-l) = NY.(kl~k~^^ (A3) 



where fc, denotes the degree of the ith node. Rewrite 



HiLi fi^i) in continuous form: 



= / Np{k)f{k)dk (A4) 



where kmin and kmax denote the minimal and maximal 
degree in RAN, respectively. Then, Equ.(A3) can be 
rewritten as: 



C 



= ,r''?^dk-2t-jf±dk (A5) 



In section 4, we have proved that p(k) = ak^"' with 7 = 3 
and a a constant, thus one can write down the expression 
that: 



C = 6a 



/•kmax 
J kmin 



^dk - 2a 



Note that l/fc3(/e - 1) = 1/(A; - 1) - l/k - l/fc^ - l/fc3^ 
one can immediately obtain the value of C as: 



/-^o/'^i"'" "i "i "i "^T ^minikmax 1)\ / \ n\ 

<^ = ^"(p— + -j— - I P (1, . _lJ 



It is clear that kmin = 3 and for sufficient large N ^ 
kmax ^ kmim and a satisfy the normalization equation: 




Therefore, a = 18 and C = f - 361n| « 0.74. 



APPENDIX B: THE AVERAGE DISTANCE 

At first, let's prove an interesting property about the 
shortest paths in RAN. Marked each node according to 
the time when the node is added to the network(see figure 
3), then we have the following lemma: 



Lemma: for any two nodes i and j, each shortest 
path from i to j{SPij) does not pass through any nodes 

k satisfying that k > maxli, j}. 
Proof. Using nodes' sequence 

1^X1^X2-^ > Xn ^ j 

to denote the shortest path from i to j of length n + 1. 
Obviously, n = is the trivial case. Suppose that n > 
and Xk = max{xi, X2, ■ • ■ ,Xn}, if VI < k < n,Xk < 
max{i,j}, then the proposition is true. 

In succession, we prove that the case Xk > max{i,j} 
would not come forth. Suppose that the triangle /^yiy2y3 
is selected when node Xk is added. Since Xk > max{i,j}, 
neither i nor j is inside the triangle Ayiy2y3. Hence the 
path from i to j passing through Xk must enter into and 
leave Ayiy2y3- We may assume that the path enter into 
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'^yi2/22/3 by node yi and leave from node then there 
exists a subpath of SPij from yi to j/2 passing through Xk , 
which is apparently longer than the direct path yi — > y2- 
Hence if SPij is the shortest path, the youngest node 
must be either i or j. ■ 

Using symbol d{i,j) to represent the distance between 
i and j, the average distance of RAN with order N, de- 
noted by L{N), is defined as: 



2a{N) 
N{N - 1) 



where the total distance is 



l<i<j<N 



(Bl) 



(B2) 



According to the lemma, newly added node will not 
affect the distance between old nodes. Hence we have: 



N 



a{N + I) ^ aiN) +J2d{i, N + 1) 



(B3) 



Assume that the node + 1 is added into the triangle 
A?/i?/22/3, then the Equ.(B3) can be rewritten as: 



N 



N 



a{N+l) = a{N)+Y^{D{i, y)+l) = a{N)+N+Y^ D{i, y) 

i=l i=l 

(B4) 

where D{i,y) = min{d{i,yi),d{i,y2),d{i,y3)}. Con- 
strict Ayiy2y3 continuously into a single node y, then 
we have D{i,y) = d{i,y). Since d{yi,y) = d{y2,y) = 
d{y3,y) = 0, the Equ.(B4) can be rewritten as: 



a{N + 1) = a{N) + A^ + ^ ^) 



(B5) 



where F = {1, 2, • • • , N} — {yi, y2, 2/3} is a node set with 
cardinality A^ — 3. 

The sum X^ier "^(^j v) be considered as the total 
distance from one node y to all the other nodes in RAN 
with order N—2. In a rough version, the sum X^ier "^l*' v) 
is approximated in terms of L{N — 2): 



J2d{i,y)^iN~3)LiN~2) 



(B6) 



Note that, the average distance L{N) increases 
monotonously with A^, it is clear that: 

(^-3)L(^-2).^^i^<^ (B7) 



n- 2 



A^ 



Combining (B5), (B6) and (B7), one can obtain the in- 
equation: 



cr(Ar + 1) < ct{N)+N + 



2a{N) 
N 



(B8) 



If (B8) is not an inequation but an equation, then the 
increasing tendency of (y{N) is determined by the equa- 
tion: 



da{N) ^ 2<j{N) 
dN N 

This equation leads to 

cf{N) = A^^lnA^ + iJ 



(B9) 



(BIO) 



where is a constant. As u{N) ^ N'^L{N), we have 
L{N) ~ InN. Which should be pay attention to, since 
(B8) is an inequation indeed, the precise increasing ten- 
dency of L may be a little slower than InA^. 
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